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Aims. We present a detailed study of the two Sun-like stars KIC 7985370 and KIC 7765 135, aimed at determining their activity level, 
spot distribution, and differential rotation. Both stars were discovered by us to be young stars and were observed by the NASA Kepler 
mission. 

Methods. The fundamental stellar parameters (v sin i, spectral type, r eff , log g, and [Fe/H]) were derived from optical spectroscopy by 
the comparison with both standard-star and synthetic spectra. The spectra of the targets allowed us also to study the chromospheric 
' j ' ■ activity from the emission in the core of hydrogen Ha and Ca n infrared triplet (IRT) lines, revealed by the subtraction of inactive 

templates. The high-precision Kepler photometric data spanning over 229 days were then fitted with a robust spot model. Model 
^0 \ selection and parameter estimation are performed in a Bayesian manner, using a Markov chain Monte Carlo method. 

Results. Both stars came out to be Sun-like (G1.5 V spectral type) with an age of about 100-200 Myr, based on their lithium content 
and kinematics. Their youth is confirmed by the high level of chromospheric activity, which is comparable to that displayed by the 
early G-type stars in the Pleiades cluster. The Balmer decrement and flux ratio of the Can-IRT lines suggest that the formation of 
the core of these lines occurs mainly in optically-thick regions that are analogous to solar plages. The spot model applied to the 
Kepler photometry requires at least seven enduring spots in the case of KIC 7985370 and nine spots in the case of KIC 7765 135 
^ ■ for a satisfactory fit of the data. The assumption of longevity of the star spots, whose area is allowed to evolve in time, is at the 

heart of our spot-modelling approach. On both stars the surface differential rotation is Sun-like, with the high-latitude spots rotating 
slower than the low-latitude ones. We found, for both stars, a rather high value of the equator-to-pole differential rotation (dQ a 0.18 
radd~') which is in contrast with the predictions of some mean-field models of differential rotation for fast-rotating stars. Our results 
are instead in agreement with previous works on solar- type stars and with other models which predict a higher latitudinal shear, 
increasing with equatorial angular velocity, that can undergo changes along the magnetic cycle. 

' Key words. Stars: activity - stars: starspots - stars: rotation - stars : chromospheres - stars: individual: KIC 7985370 and KIC 7765135 

' - X-rays: stars 
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f^*) . 1 . Introduction belts, absence or a different behaviour of activity cycles, highly 
' energetic flares, etc. These differences are likely related to the 
In the Sun, magnetic activity is thought to be produced by a dynamo mec hanism, which is operating in rather different con- 
global-scale dynamo action arising from the coupling of con- ditions in young stars mainly regar ding their rotation rate and 
vection and rotation ( [Parker] 1955HSteenbeck et al. | 1966|) . Young internal stmcture . Understanding the properties of young suns 
V ; Sun-like stars are rotating faster than the Sun and display a much and; p articu i ar i y , their activity and rota tion, is crucial to trace 
^ . higher level of magnetic activity at all atmospheric layers, which the Sun and its environment back to the first evolutionary stages. 
_ CS is likely due to a stronger dynamo action. They also show differ- In fact the properties of the magnetoconvection in these stars 
ent manifestations of activity compared to the Sun, such as big- seems tQ be strongly innue nced by the Q -effect that pr oduces 
ger and long-living spots in their photospheres, active longitude characteristic "wreaths" on large scales (iNelson et al.l l20ll . 

Although strong latitudinal variations of the differential rota- 
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e-mail: hefroehlich@aip.de wind balance and geostrophy, the results of numerical simula- 

* Based on public Kepler data, on observations made with the Italian tions seem to be strong i y dependent on the Reynolds number of 

Telescopic Nazionale Galileo (TNG) operated by the Fundacion Galileo ^ flow The situation is at variance wim the dynamo action in 

Galilei - INAF at the Observatono del Roque del los Muchachos, La • , . u , . , 

, ■ ii j i-i-, i main sequence solar-type stars, where the role or the tachochne 

Palma (Canary Islands), on observations collected at the 2.2-m tele- . . , . , . Jr . , „ nr-r . „ 1 

scope of the Centro Astronomico Hispano Aleman (CAHA) at Calar ismstead essential in producing the a effect (|Dikpati & Gilman| 

Alto (Almeria, Spain), operated jointly by the Max-Planck-Institut fur I2Q0J1 | Bonanno et al. || 2002|; | Bonanno || 2012D . 

Astronomie and the Instituto de Astrofisica de Andalucfa (CSIC), and It is still unclear whether a strong latit udinal differential rota- 

on observations collected at the Catania Astrophysical Observatory tion is common among fast-rotating stars. iMarsden et all (1201 ll) 

(Italy). report values of the absolute differential rotation dQ in the range 
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0.08-0.45 radcr 1 for a sample of rapid rotators similar to and 
slightly more massive than the Sun. Despite the spread of val- 
ues, it seems that d£l is in any case larger than in the Sun. The 
measures of absolute differential rotation in a large sample of F- 
and early G-type stars through the Fouri er transform technique 
dReiners & Schmittll2003l |R einers I I2006T) show no indication of 
the decrease in this parameter with the rotation period, rather 
the highest values of dQ are encountered for periods between 2 
and 3 days. On the other hand, some recent calculations predict 
a moderate differential rotation, comparable to that of the Sun, 
also f or a Sun-like star rotating 20 times faster (e.g.. lKiiker et al.l 
l20ll . 

Moreover, d£l seems also to be a function of the stellar mass 
for main-sequence stars, increasing with the ir effective temper- 
ature, as shown, e.g., bv lBarnes et al. 1 (12001 . One of the largest 
values of differential rotat ion for a star noti ceably cooler than 
the Sun was found by us dFrasca et al.ll2.01 ll hereafter Paper I) 
for KIC 8429280, a 50Myr-old K2-type star, from the analysis 
of the light curve collected by the NASA Kepler spacecraft. 

T he highly precis e photometry of Kepler dBorucki et al.l 

l2010tlKoch et alj2010h coupled with the long and virtually unin- 
terrupted coverage makes these data unique for the study of pho- 
tospheric activity and differential rotation in late-type stars, as 
we have shown in our first work based on the analysis of Kepler 
data of spotted stars (Paper I). 

However, whether star spots are indeed the best tracers of the 
surface rotation o r not is still a matter of deb ate (for a different 
point of view, see lKorhonen & Elstnerll20Tlh . 

As for KIC 8429280 (Paper I), the two new targets 
KIC 7985370 (HD 189210 = 2MASS J19565974+4345083 = 
TYC 3149-1571-1) and KIC 7765135 (2MASS J19425057 
+4324486 = TYC 3148-2163-1) were selected as active stars 
from their optical variability and from th e cross-corr e lation 
of th e ROSAT All-Sky Survey (RASS; IVoges et al.l [l999l 
200C ) with Tycho and Hipparcos catalogues (iPerrvman et al.l 
19971) . With V = lOT'O and 11™8, respectively, both stars 



are relatively bright ones in the Kepler field o f view. Both 
of the m were recently reported as variable by iPigulski et al.l 
(I2009h . who searched for bright variable stars in the Kepler 
field of view with ASAS3-North station. The variability of 
K IC 7985370 could be due to a rotational modulation according 
to lUvtterhoeven et all d201 ll) that rely on the first two quarters 
of Kepler data. The Kepler light curves readily show these stars 
as rotationally variable with a period of about 2-3 days, which 
is typical of G-type stars in the P leiades cluster (age « 130Myr, 
iBarrado v Navascues et al.ll2004l) . The estimates of their atmo- 
spheric parameters reported in the Kepler Input Catalog (KIC), 
which ar e based on Sloan photome try (for a revised temperature 
scale cf. iPinsonneault et al] 1201 2|) . suggested to us that these 
objects were similar to the Sun. 

The analysis of the optical spectra collected by us confirmed 
that the stars are nearly identical to the Sun, but much younger 
and as such deserving a detailed investigation. Applying the 
same techniques as in Paper I, we determined the basic stel- 
lar parameters (Sect. 13. lV the chromospheric activity (Sect. l3~2l ). 
and elements abundances (Sect. 13.31 ). The kinematics of these 
stars is shortly discussed in Sect. 13. 51 The Bayesian approach 
to spot modelling applied to the Kepler light curves is described 
in Sect.|4]and the results are presented and discussed in Sect. [5] 
with particular relevance to the issue of differential rotation. 



2. Ground-based observations and data reduction 

2.1. Spectroscopy 

Two spectra of KIC 7985370 with a signal-to-noise ratio S/N 
of about 60-70 were collected at the M. G. Fracastoro station 
(Serra La Nave, Mt. Etna, 1750 m a.s.l.) of the Osservatorio 
Astrofisico di Catania (OAC, Italy) in July and October 2009. 
The 91 -cm telescope of the OAC was equipped with FRESCO, a 
fiber-fed echelle spectrograph that, with the 300-lines/mm cross- 
disperser, covers the spectral range 4300-6800 A with a resolu- 
tion = A/AA * 21000. 

Another spectrum of KIC 7985370 with S/N* 120 was taken 
with SARG, the echelle spectrograph at the Italian Telescopio 
Nazionale Galileo (TNG, La Palma, Spain) on 2009 August 12 
with the red grism and the slit width of 078. This spectrum, 
covering the 5500-1 1 000 A wavelength range, has a resolution 
R * 57 000. 

The spectrum of KIC 7765 135 was taken on 2009 October 
3 with the Fiber Optics Cassegrain Echelle Spectrograph 
(FOCES: IPfeiffer et al.ll 19981) at the 2.2-m telescope of the Calar 
Alto Astronomical Observatory (CAHA, Almeria, Spain). The 
2048x2048 CCD detector Site#ld (pixel size = 24 /mi) and the 
slit width of 400 /mi give rise to a resolution R ss 28 000 and a 
S/N* 100 in the red wavelength region with a 30-minutes expo- 
sure time. 

Spectra of radial (RV) and rotational velocity (v sin i) stan- 
dard stars (Table[TJ, as well as bias, flat-field, and arc-lamp ex- 
posures were acquired during each observing run and were used 
for the data reduction and analysis. 

The data reduction was performed with the Echelle task of 
the IR AF3 package (see, e.g jFrasca et al.l2010tlCatanzaro et al.l 
I2010L for details). 



Table 1. Radial/rotational velocity standard stars. 



Name 


Sp. Type 


RV 


v sin f 


Notes'* 






(kms" 1 ) 


(kms" 1 ) 




HD 10307 


G1.5V 


4.0 6 * 


2.1 


v sin i 


HD 10700 


G8V 


-17.1 fc 


0.9 


v sin i 


HD 157214 


GOV 


-19.2 b 


1.6 


v sin i 


HD 182572 


G8 IV 


-100.35" 


1.9 


RV,v sin i 


HD 187691 


F8V 


0.0" 


2.8 


RV 



N otes. M IUdrvetal.1 d!999h . w iNordstrom etaD d2004l) . 
(c) lGlebocki & Gnacinskil d2005l) . ( ~j?V = radial velocity standard star; 
v sin i = standard for rotational velocity. (,) Variable? 



2.2. Photometry 

The focal-reducer CCD camera at the 91 -cm telescope of OAC 
was used, on 2009 December 10, to perform standard photome- 
try in the Johnson-Cousins B, V, Rc, and Iq bands. The data were 
reduced following standard steps of overscan region subtraction, 
master-bias subtraction, and division by average twilight flat- 
field images. The BVRqIc magnitudes were extracted from the 
corrected images through aperture photometry performed with 



1 IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of the Universities 
for Research in Astronomy, inc. (AURA) under cooperative agreement 
with the National Science Foundation. 
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Table 2. Stellar parameters. The data of the upper part is from 
the literature. 







KIC 7985370 


KIC 7765135 


RA 


(J2000) 


19 h 56 


m 59!74 


19 h 42 


"5^58 


DEC 


(J2000) 


+43°45'08'.'3 


+43°24'48'.'7 


X-ray count" 


[ct/s] 


0.033 


± 0.008 






J" 




8.789 


± 0.020 


10.414 


± 0.025 


H» 




8.507 


+ 0.063 


10.113 


+ 0.020 


K s b 




8.419 


± 0.033 


10.028 


± 0.014 


B 




10.59 


± 0.09 




-t- n oq 

± U.U7 


V 




9.98 


± 0.08 


1 1 89 


It u.uo 


R c 




9.57 ± 0.06 


1 1 98 

1 1 .La 


± u.uo 


Ic 




9.33 


± 0.06 


1 Q7 


It u.uu 


Sp. Type 




G1.5V 


v 1 1 


S V 

^' V 


T<sB 


[K] 


5815 


+ 95 


5835 


± 95 






4.24 


+ 0.12 


4.34 


± 0.12 


[Fe/H] 




-0.05 


± 0.10 


0.04 


± 0.05 


v sin i c 


[km s ] 


18.2 


± 1.3 


21.4 


± 1.1 


v sin i d 


[km s ] 


17.9 


± 1.3 


21.9 


± 1.3 




[Kill S J 


-24.0 


± 0.3 


-20.0 


+ 0.2 


W(U) 


[mA] 


155 


± 20 


160 


± 20 


log N(uy 




2.87 


± 0.10 


2.93 


± 0.10 


Age 


[Myr] 


100 


- 200 


100 


- 200 


Mass 


[M Q ] 


1.15 


± 0.10 


1.15 


± 0.10 


Distance 


[pc] 


113 


± 15 


245 


+ 30 


V 


[km s- 1 ] 


-8.63 


+ 0.86 


-5.61 


± 1.73 


V 


[km s- 1 ] 


-22.03 


± 0.37 


-18.14 


± 0.57 


w 


[km s" 1 ] 


-6.74 


± 0.84 


-7.86 


± 1.93 



Notes. (n) From t he ROSAT All-Sky Survey Faint Source Catalogue 
l Voges et al ll2000h . ( *> JHK S magnitudes are from 2MASS catalogue 
i Cutri et al.ll2003h . (c) From ROTFIT code. w From FWHM-v sin i cali- 
bration. (e) Calculated by means of a spectral synthesis based on ATLAS9 
atmospheric models (Sect. l3.3V 



DAOPHOT by using the IDL0 routine APER. Standa rd stars 
in the cluster NGC 7790 and in the field of BL Lac (IStetsonl 
|2000T) were used to calculate the zero points and the transfor- 
mation coefficients to the Johnson-Cousins system. Photometric 
data derived in th is work and JHK S magnitudes from 2MASS 
dCutri et al.ll2003l) are summarized in Table |2] 

3. Target characterization from ground-based 
observations 

3.1. Astrophysical parameters 

The analysis of the high-resolution spectra was aimed at measur- 
ing the radial (RV) and projected rotational velocities (vsin/), 
performing the MK spectral classification, and deriving ba- 
sic stellar parameters like effective temperature (T e g), gravity 
(logg), and metallicity ([Fe/H]). 

We have computed the cross-correlation functions (CCFs) 
with the IRAF task fxcor adopting spectra of late-F and G- 
type slowly-rotating stars (Table[T), which were acquired with 
the same setups and during the same observing nights as our 
targets. The latter were used both for the measure of RV and 
the determination of vsin; through a calibration of the full- 
width at half maximum (FWHM) of the CCF peak as a func- 
tion of the v sin of t heir artificially broadened spe ctra (see, e.g., 
iGuillout et aill2009l: iMartfnez-Arnaiz et al.ll2010l) . We averaged 

2 IDL (Interactive Data Language) is a registered trademark of ITT 
Visual Information Solutions. 



the results from in dividual echelle orders as described, e.g., in 
iFrasca et alJ (I2010I) to get the final RV and v sin i of our young 
Sun-like stars. 

For KIC 7985370, we measured, within the errors, the same 
RV in the FRESCO spectra taken on 2009 July 10 (RV = -23.9+ 
0.4 kms" 1 ) and on 2009 October 4 (RV = -24.2 + 0.4 kms" 1 ), 
while a value of RV = -21 .8 + 0.2 km s _1 has been derived from 
the SARG spectrum of 2009 August 12. This could indicate a 
true RV variation just above the 3-cr confidence level. However, 
much more spectra are needed to prove or disprove such a RV 
variation and to make sure that the star is a single-lined spectro- 
scopic binary. 

For KIC 7765 135, we only acquired one spectrum from 
which we derived an RV of -20.0+0.2 km s . 

Since we found no indication for binarity (visual or spectro- 
scopic) in the literature for both stars, we consider them as being 
single stars or, at most, single-lined spectroscopic binaries for 
which the eventual companions are too faint to give a significant 
contribution neither to the optical spectrum nor to the Kepler 
photometry. 

For the SARG spe ctrum, we used the ROTFIT code 
(IFrasca et alj|2003l 120061) to evaluate T eS , logg, [Fe/H], and re- 
determine v sin i, adopting a library of ELODIE Archive spectra 
of standard stars, as described in Paper I. To match the lower res- 
olution (Afoces = 28 000) of the FOCES spectrum, the ELODIE 
templates (/?elodie = 42 000) were convolved with a Gaussian 

kernel of FWHM = Ajl/Rp QCES - before running 

ROTFIT. We applied the ROTFIT code only to the echelle orders 
with a fairly good S/N, which span the ranges 5600-6700 A and 
4300-6800 A for SARG and FOCES spectra, respectively. The 
spectral regions heavily affected by telluric lines were excluded 
from the analysis. The adopted estimates for the stellar parame- 
ters and the associated uncertainties come from a weighted mean 
of the values derived for all the individual orders, as described 
in Paper I. We also applied this procedure to the FRESCO spec- 
tra of KIC 7985370, obtaining stellar parameters in very close 
agreement with those derived from the SARG spectrum. 

ROTFIT also allowed us to measure v sin i by matching to 
the observed spectrum the spectra of slowly-rotating standard 
stars artificially broadened at an increasing vsin; and finding 
the minimum of x 1 - F° r this purpose, we used the standard stars 
listed in TableQ] because their spectra were acquired with the 
same instrumental setup as that used for our targets, avoiding 
the introduction of any systematic error caused by a different 
resolution. 

We obtained nearly the same values of the astrophysical pa- 
rameters for both targets. They turn out to be very similar to the 
Sun, both in spectral type and effective temperature, but consid- 
erably younger, as testified by the high lithium content in their 
photospheres (Sect. l3.3l l. Their average values and standard er- 
rors are reported in Table 

Based on these parameters and the age estimate (Sects. l3~3l 
and 13.5b . compatible with zero-age main-sequence (ZAMS) 
stars, we can roughly estimate their mass and radius as M — 
1.15 ± O.lOMfp and^ = 1.1 + 0.1 Rq from evolutionary tracks 
(e.g. JSiess etai]l2000l) . 

3.2. Chromospheric activity 

The high level of magnetic activity in the chromosphere is dis- 
played by the Ca n IRT and Ha lines that are strongly filled-in by 
emission (Figs.[TJand[2]i. The Hei D3 (/15876A) absorption line 
is clearly visible in the spectra of our two young suns (Figs.[TJand 
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Hel Ha Call AA8498.8542 Call A8662 





5860 5865 5870 5875 5880 5885 
A (A) 



6550 



6560 6570 
X (A) 



6580 



8500 8520 8540 
X (A) 



8560 



8650 8660 8670 8680 
X (A) 



Fig. 1. Top of each panel: Observed, continuum-normalized SARG spectrum of KIC 7985370 (solid line) in the Hei D3, Ha, and 
Can IRT regions together with the non-active stellar template (dotted red line). Bottom of each panel: Difference between observed 
and template spectra. The residual Ha profile is plotted shifted downwards by 0.2 for the sake of clarity. The hatched areas represent 
the excess emissions (absorption for He 1) that have been integrated to get the net equivalent widths. 




Call AA8498.8542 



Call A8662 




5860 5865 5870 5875 5880 5885 
X (A) 



6550 6560 6570 6580 
X (A) 



Fig. 2. As Fig.Q] but for the FOCES spectrum of KIC7765135. 



8500 8520 8540 8560 
X (A) 



8650 8660 8670 8680 
X (A) 



0. This implies the presence of an upper chromosphere because 
a temperature of at least 10 000 K is required for its formation. 

We estimated the chromospheric emiss ion level with the 
"spectral subtraction" tec hnique (see, e.g., iFrasca & Catalanol 
Il994t iMontes et al.lll995t) . This procedure is based on the sub- 
traction of a reference "non-active" template made with ob- 
served spectra of slowly-rotating stars of the same spectral type 
as the investigated active stars, but with a negligible level of 
chromospheric activity and an undetectable lithium line. The 
equivalent width of the lithium and helium absorption lines (Wu 
and Wu e , respectively) and the net equivalent width (W em ) of the 
Can IRT and Ha lines were measured from the spectrum ob- 
tained after subtracting the non-active template from the target 
and then by integrating the residual emission (absorption for the 
lithium and helium lines) profile, which is represented by the 
hatched areas in Figs.[T]and[2] 

We noted a significant variation of the Ha emission in the 
three spectra of KIC 7985370, with ranging from 160 to 

264 mA (TableE). 

The radiative losses in the chromospheric lines were eval- 
uated by multiplying the average W em by the continuum sur- 
face flux at the wa velength of the line, in the same way as in 
IFrasca et al.l (120 1 0b . The net equivalent widths and the chromo- 
spheric fluxes are reported in Table [3] 

The HB net equivalent width could be measured only for 
KIC 7765 1 35 because the SARG spectrum of KIC 7985370 does 



not include that wavelength region. In the FRESCO spectra, it is 
not possible to apply the subtraction technique, due to the very 
low S/N ratio in that spectral region, which prevented us from 
detecting such a tiny filling of the line. 

For KIC 7765 135 we measured a Balmer decrement 
Ffta/Fiip = 4.3 + 2.5 that is slightly larger than the values 
in the range 2-3 found by us for K IC 8429280 (Paper I) and 
for HD 171488 (IFrasca et al.ll2010l) . Accounting for the error, 
this ratio is however still compatible with optically-thick emis- 
sion b y atmospheric features simil ar to solar and stellar plages 
(e.g., iBuzasilfl989T: IChesterifl99ll) rather than to prominence - 
like structures, for which a Balmer decrement of the order of 10 
is exp ected (e.g., lLandman & Mongillolll979l : lHall & Ramsevl 
[199~ 



The Can-IRT flux ratio, ^8542/^8498 = 1.4 ± 0.9, is in- 
dicative of h igh optical dep ths and is in the range of the val- 
ues found by IChesterl d!991l) in solar plages. The optically-thin 
emission from solar prominences gives rise instead to values of 
-F8542/-F8498 ~ 9. We measured a value of / r 8542/f 8498 = 1-6+1.0 
for KIC 7985370, that is almost the same as KIC 7765 135. 

The chromospheric emission in both stars appears thus to be 
mainly due to surface regions similar to solar plages, while the 
quiet chromosphere and eventual prominences play a marginal 
role. 

The subtraction technique also allowed us to measure the 
lithium equivalent width cleaned up from the contamination 
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Table 3. Line equivalent widths and chromospheric fluxes. 



Line 


Date 


UT mid 




Error 


Flux 




(yyyy/mm/dd) 


(hh:mm) 


(niA) 


(mA) 


(erg cm -2 s _1 ) 






KIC 7985370 








2009/07/10 


22:29 


160 


40 


1.23xl0 6 




2009/08/12 


04:18 


264 


30 


2.02xl0 6 


" 


2009/10/03 


19:17 


202 


35 


1.55xl0 6 


He i D 3 


2009/08/12 


04:18 


-50" 


20 




Can ,18498 


» .. 




263 


55 


1.33x10 s 


Can/18542 


» » 




439 


55 


2.16x10 s 


Can /18662 


" " 




337 


60 


1.66x10 s 






KIC 7765 135 






Ha 


2009/10/03 


20:15 


326 


60 


2.48x10 s 


W 






54 


30 


0.58x10 s 


Hei D 3 






-60" 


25 




Caii/18498 






258 


45 


1.28x10 s 


Can /18542 






366 


45 


1.76x10 s 


Can /18662 






328 


65 


1.58x10 s 



Notes. <n> The minus sign for the VK cm indicates a residual absorption 
(He i line). 



of Fei A 6707.4 A line, which is strongly blended with the 
nearby Li i line (Fig. [3]). The lithium equivalent width measured 
for KIC 7985370 (Wu = 155 + 20 mA) and for KIC 7765 135 
(Wu - 160 ± 20 mA), jus t above the Pleiades upper enve- 
lope dSoderblom et al.l 1993 a). translates into a very high lithium 
abundance, log N( Li) = 3.0 ± 0.1, for both star s, as results from 
the calibrations of lPavlenko & Magazzul(ll996l) . The determina- 
tion of lithium abundance was refined in Sect. l3.3l bv me ans of a 
spectr al synthesis based on ATLAS9 atmospheric models dKuruczl 
119931) . 

3.3. Abundance analysis 

Since the spectral features in our two stars are relatively broad, it 
is difficult to find unblended lines for measuring their equivalent 
widths. To overcome this problem, the photospheric abundances 
were thus estimated by matching a rotationally-broadened syn- 
thetic spectrum to the observed one. For this purpose, we divided 
each spectrum into 50 A-wide segments, which were analyzed 
separately. The spectral ranges 4450-8650 A and 5600-7300 A 
were used for the spectra of KIC 7765 135 and KIC 7985370, 
respectiv ely. The synthetic line profiles were c omputed with 
SYNTHE (lKurucz&Avrettlll98ll) using ATLAS 9 (lKuruczHl993h 
atmospheric models. All models were calculated using the so- 
lar opacity distribution function and a microturbulence velocity 
£ = 2 km s _1 . For each segment, the abundance was determined 
by^ 2 minimiz ation. We used the spectral line list and atomic pa- 
rameters from ICastelli & Hubrig I (12004b that are updated from 
lKurucz&Beim995h . 

The error in the abundance of one particular element was 
taken as the standard deviation of the mean of the abundances 
calculated for each segment. For elements whose lines occurred 
in only one or two segments, the error in the abundance was 
evaluated by varying the effective temperature and gravity within 
their uncertainties (as given in Table [2) and computing the abun- 
dance for r e ff and logg values within these ranges. We found 
a variation of » 0. 1 dex in abundance due to temperature varia- 
tion, but no significant variation when log g was varied. Thus the 
uncertainty in temperature is probably the main source of error 
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Fig. 3. Observed spectra (dots) of KIC 7985370 and 
KIC7765135 (shifted downwards by 0.4) in the Lii A 6707.8 A 
region together with the synthetic spectra (full lines). 

in the abundance estimation. In Fig. [3] we show an example of 
the matching between synthetic and observed spectra. To deter- 
mine the Li abundance, we used th e Li i AblOl A line, t aking 
into account the hyperfine structure (lAndersen et al.|[T984b . The 
abundances are given in Table [4] 

Fig. [4] shows a comparison between the abundances derive d 
for our targets and the solar values from iGrevesse et al 1 (120101) . 
With the exception of lithium, all the other abundances are rather 
similar to those measured for the solar photosphere. In particular, 
we note that KIC 7765 135 is a:0.3 dex more metallic than the 
Sun, which confirms the results of ROTFIT (Tabled. 

As far as the lithium is concerned, we converted the 
abundance from log N/N tot to logN/N{H) on a scale where 
logJV(H) = 12, finding logJV(Li) = 2.93 ± 0.10 and 2.87 ± 
0.10 for KIC7765135 and KIC7985370, respectively. These 
values are smaller than, but marginally consistent with, those 
obtained from the equivale nt widths and the calibrations of 
iPavlenko & Magazzul d 19961) and reported in Sect.l3~S1 

For early G-type stars, like our targets, the lithium abun- 
dance does not give a very strong constraint on their age. 
Indeed, their log N(Li) remains close to the initial lithium abun- 
dance for Population I stars (log /V(Li)=3. 1-3.3) during their 
early life. It has been shown that stars with temperatures in 
the range T e g = 5900 ± 150 K display a significant lithium 
depl etion (log/V(Li) = 2.9 ) only after about 150 Myr (see, 
e.g., ISestito & Randichll2005l) . Thus, we can estimate an age of 
the order of 100-200 Myr for both stars with a lower limit of 
roughly 30-50 Myr, being the lithium abundance below the up- 
per envelope of the or Per (age » 50 Myr) and IC 2602 (age » 
30 Myr) clusters dMontes et al .11200 1 bl: [Sestito & Randichl2005l) . 
Moreover, the absence of a strong near-IR excess in their spectral 
energy distributions (see Sect. 13. 4t and of accretion signatures in 
the spectrum, allows us to exclude an age of a few million years 
for both stars that should already be in a post-T Tauri phase (age 
> 10-20Myr). 
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Table 4. Abundances inferred for our two stars expressed in the 
form log N/Nto t . 
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Li 


-9.10 ± 


0.10 
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C 


-3.67 ± 


0.10 


-3.29 ±0.14 


O 






-3.10 ±0.10 


Na 


-5.77 ± 


0.16 


-5.68 ±0.12 


Mg 


-4.67 ± 


0.11 


-4.77 ±0.13 


Al 


-5.48 ± 


0.17 
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Si 
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Ca 
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Ti 
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-6.82 ± 0.09 
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Fig. 4. Abundances found for KIC 7765 135 (open circles) and 
for KIC 7985370 (filled circles) related to the solar ones. 



3.4. Spectral energy distributions 

The spectral energy distributions (SEDs) from the optical to the 
near-infrared (IR) domain were obtained by merging our stan- 
dard UBVRqIc photometry (Table [2]i with JHK S magnitudes 
from the 2MASS catalogue dCutri et al.ll2003l) . 

We used the grid of NextGen low-resolution synthetic 
spectra, with log;? - 4.0 and 4.5 and solar metallicity by 
lHauschildt eTaT] d 19991) . to perform a fit to the SEDs. The ef- 
fective temperature (T e s) was kept fixed to the value derived 
with ROTFIT (Table [2]i for each target. The interstellar extinc- 
tion (Ay) was evaluated from the distance according to the rate 
of 0.8 mag/kpc found by iMikolajewska & Mikolaiewskil Jl98dh 
for the sky region aro und CI Cyg (very close to our targets). The 
ICardelli et alld 19891) extinction law with Ry — 3.1 was used for 
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Fig. 5. Spectral energy distributions (dots) for KIC 7985370 (top 
panel) and KIC 7765135 (bottom panel). The NextGen synthetic 
spectrum at r e ff = 5800 K scaled to the star distance is over- 
plotted with continuous lines in each box. 



evaluating the extinction in the other bands. Finally, the angular 
diameter (<p) that scales the synthetic surface flux over the stellar 
flux at Earth, was left free to vary. The best solution was found 
by minimizing the^- 2 of the fit to the BVRqIcJ data, which are 
dominated by the photospheric flux of the star and are normally 
not significantly affected by infrared excesses. The angular di- 
ameter derived for KIC 7985370, <p = 0.0895 mas, implies a dis- 
tance d = 113 ± 15 pc, if we adopt a ZAMS radius of 1.1 Rq 
(SectEQ. For KIC 7765135 we found = 0.0413 mas, which 
corresponds to d — 245 + 30 pc, adopting again 1 . 1 Rq. The error 
in the distance is estimated taking into account both a 0. I-Rq un- 
certainty in the stellar radius and the temperature error of about 
100 K (see Table©. 

As apparent in Fig. [5] the SEDs are well reproduced by the 
synthetic spectrum till the K s band both for KIC 7985370 and 
KIC 7765 135 and no excess is visible at near-IR wavelengths. 
This strengthens the spectroscopic determination of the effective 
temperature and ensures that these stars have by far gone over 
the T Tauri phase, during which a thick and dense accretion disk 
gives rise to a conspicuous infrared excess. The absence of mid- 
and far-infrared data does not allow us to exclude the presence of 
thinner "d ebris" disks, like those found in a few Pl eiades solar- 
type stars dStauffer et aLlfeOOStlGorlova et al.ll20Q6b . 
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Fig. 6. (U, V) and (W, V) planes of our two young suns. We plotted the average position of each young stellar kinematic group with 
different symbols. The dotte d line in the left pa nel demarcates the locus of the young-disk population (age < 2 Gyr) in the solar 
neighborhood, as defined bv lE^n1(fT98l[T98l . 



3.5. Kinematics 

We used the radial velocities determined in Sect. 13.11 and the 
distances estimated in Sect. 13.41 together with measu rements of 
prope r motions taken from the Tycho-2 catalogue (H0g et al. 
l2000h . to calculate Galactic space-veloci ty components (U, V , 
W), following the procedures described in lMontes et al .1(1200 lab . 
The values, listed in Table [2] are given in a right-handed coor- 
dinate system (positive toward the Galactic anti-center, in the 
Galactic rotation direction, and toward the north Galactic pole, 
respectively). In the {U, V) and (W, V) Bottlinger diagrams 
(Fig. [6}, we have plotted the locus of our targets. Different 
sym bols represent the ce ntral position given in the literature 
(see iMontes et alj 1200 lab of the five youngest and best docu- 
mented moving groups (MGs) and superclusters (SCs): namely 
the Local Association (LA) or Pleiades MG (20-150 Myr), the 
IC2391 SC (35-55 Myr), the Castor MG (200 Myr), the Ursa 
Major (UMa) MG or Sirius supercluster (300 Myr), and the 
Hyades SC (600 Myr). 

The space velocities of these two Sun-like stars are consis- 
tent with those of the young-disk population (Fig. |6). Based 
on t wo different st atistic al methods (for further information, 
see IKlutsch et all l[20Tot) and the paper in preparation by 
Li" 



IKlutsch et alJ d2012h . we determined their membership proba- 
bility to each of the five aforementioned young stellar kine- 
matic groups. KIC 7985370 and KIC 7765135 fall in the Local 
Association (LA) locus. With a probability of more than 70 %, 
they turn out to be highly likely members of this moving group. 
Furthermore, such a result is absolutely in agreement with the 
age derived from the lithium abundance (Sect. 13.3b . 



4. Spot modelling of the Kepler light curves 

4.1. Photometric data 

All the available public Ke pler long-cadence time-series (At » 
30 min. Uenkins et al.ll2010l) . spanning from 2009 May 2 to 2009 
December 16, was analyzed. It covers altogether 229 days and 



corresponds to the observing quarters 0-3 (Q0-Q3), with the 
largest gap, about 4.5 days, appearing between Ql and Q2. 

To remove systematic trends in the Kepler light curves asso- 
ciated with the spacecraft, detector, and environment, and to pre- 
pare them for the analysis of star spots that we will describe be- 
low, we used the software kepcotreni0. This procedure is based 
on Cotrending Basis Vectors (CBV), which are calculated (and 
ranked) through singular value decomposition and describe the 
systematic trends present in the ensemble flux data for each CCD 
channel. We used the first two basis vectors for Q0 data, while 
from three to five CBV were adopted for the correction of longer 
data sets such as Ql, Q2, and Q3. 

In order to check how the data rectification accomplished 
with kepcotrend is reflected in the outcome, the spot modelling 
has been done twice: with the rectified data (Case A) as well as 
with the original data (Case B). 

The power spectra of the K epler time-ser i es, cl eaned by 
the spectral window according to iRoberts et al.l (119871) . are dis- 
played in Fig. UJ The lower panel of Fig. [7] clearly shows 
two main peaks for KIC 7765135, which are close in fre- 
quency (0.391 and 0.414 d~'). The corresponding periods are 
2.560±0.015 and 2.407+0.014 days, respectively. The period 
errors are from the FWHM of the spectral window. The low- 
amplitude peaks at frequency of «0.8 d -1 are overtones of the 
two main peaks. As visible from the upper panel of Fig. [7j the 
structure of the peaks for KIC 7985370 is more complex, with 
the maximum corresponding to 2.856±0.019 days and a second 
peak, blended with the first one on its low-frequency side, at 
2.944 days. A third small peak corresponding to P — 3.090 days 
is also visible. 

Such a double- o r multiple-peaked p eriodogram hints at dif- 
ferential rotation. As lLanza et aD (1 19941) predicted, a photomet- 
ric time series, if sufficiently accurate (AF/F = 10~ 5 - 10~ 6 ), 
may reveal a Sun-like latitudinal differential rotation. 

An estimate of the inclination of the rotation axis with re- 
spect to the line of sight is very useful to constrain the spot 



3 http://keplergo.arc.nasa.gov/ 
ContributedSof twareKepcotrend . shtml 
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Fig. 7. Cleaned periodograms of the Kepler Q0+Q1+Q2+Q3 
time series for KIC 7985370 (upper panel) and KIC 7765 135 
(lower panel). 



model. With vsinz, stellar radius R, and rotation period P 
known, the inclination of the rotation axis follows from 



sin; 



(v sin i) ■ P 
2^R 



(1) 



In the absence of an accurate parallax value the stellar radius 
cannot be derived from the effective temperature and luminos- 
ity. If we adopt the radius for a ZAMS star with the effective 
temperature of our targets (T e g = 5800 K), R as 1.1 Rq, we get 
sin i = 0.967 (i = 75° ± 15°) for both stars. However, as stated in 
Sect. 13.31 the lithium content cannot provide a firm lower limit 
for the ages of these stars, which could also be as young as a few 
10 Myr (post-T Tauri phase). Thus, allowing for such a young 
age, a Teg = 5800 K is reached by a star of 1.5 Mq at 10 Myr 
wi th a radius of abou t 2 Rq according to the evolutionary tracks 
bv lSiess et al ] (120001) . In this case, an inclination of about 30° is 
deduced. 



4.2. Bayesian photometric imaging 

The method is basically that of Paper I. However, as the light 
curves are now significantly longer, the introduction of further 
free parameters was inevitable. The latitudinal dependence of 
rotation frequency Q(B) (Eq.[2ji now contains a sin ;8-term and, 
more important, the prescription for spot area evolution is much 
more detailed. Furthermore, the likelihood function (Eq.|3j is 
generalized by taking into account an unknown linear trend in 
the data. To tackle the problem of strong correlations between 
some parameters, an essential new ingredient is the usage of an 
orthogonalized parameter space, where the steered random walk 
of the Markov chains is performed. For the sake of clarity, due to 
these new features, the basics of the method are explained in this 
section, although they can also found in Paper 1 . A ful l account 
of the method must be deferred to a forthcoming paper dFrohlichl 
l20l2h . 

A light-curve fitting that represents spots as dark and circular 
regions has the advantage of reducing the dimensionality of the 
problem and to promptly provide us with average parameters 
(area, flux contrast, position, etc.) for each photospheric active 
region. Of course, there are other techniques, which are based 



on different assumptions, to reconstruct surface features photo- 
metrically as the in version of Kepler light curves done by, e. g., 
iBrown et all d201 ll) . 

Our aim is to present a low-dimensional spot model, with 
few spots only, that fits reasonably well the data regardless to 
very low-amplitude details that require a high degree of com- 
plexity. In a Bayesian context this claim could be even quan- 
tified. Ideally, one should estimate the so-called evidence, the 
integral over the posterior probability distribution. It would pro- 
vide a measure of the probability of a n-spot model and, there- 
fore, allow one to constrain the number of spots n that are re- 
ally needed. For numerical reasons we are compelled to resort 
instead t o the less deman ding Bayesian Information Criterion 
(BIC) bv lSchwarzl dl978h . The latter or any other related crite- 
rion expresses Occam's razor in mathematical terms without the 
need to compute the evidence. Unfortunately, we have to admit 
that - due to the unprecedented accuracy of the Kepler data - we 
could not reach this goal with only seven or nine spots. There is 
obviously more information in the data than our most elaborate 
model is able to account for. 

Dorren's (1987) analytical star-spot model, generalized to a 
quadratic limb-darkenin g law, was used. The two co efficients are 
taken from the tables of Claret & Bl oemen l d2oTTb for a micro - 
turbulence velocity of £ = 2kms _1 and are used for both the 
unperturbed photosphere and the spots. 

Four parameters describe the star as a whole: One is the co- 
sine of the inclination angle i. Three parameters (A, B and C) de- 
scribe the latitudinal dependence of the angular velocity. With B 
being the latitude value, the angular velocity Q is parameterized 
by a series expansion using Legendre polynomials: 

CIQ3) = A + 3B{5 sin 2 /? - l) /2 + 

+ c(315sin 4 /3-210sinV+ 15)/8. (2) 

The equatorial angular velocity is £2 e q - A - 3B/2 + \5C/% and 
the equator-to-pole differential rotation AQ. = 15,6/2 + 105C/8. 
In the case of the equator rotating faster than the poles AQ. is 
negative. In what follows the minus sign is suppressed, and only 
the absolute value |dQ| is given. Both stars are definitely rotating 
like the Sun. 

As in Paper 1, all star spots have the same intensity k rela- 
tive to the unspotted photosphere and are characterized by two 
position coordinates (latitude and initial longitude) and by their 
radius. All these are free parameters in the model. Other spot pa- 
rameters are the rotation period, which defines the spot longitude 
at any time and is tied to the latitude via Eq.|2] The hemisphere, 
to which a spot belongs to, is to be found by trial and error. 
Further parameters describe the spot area evolution. 

As our photometric analysis mainly aims at estimating the 
level of surface differential rotation, our focus is on long-lived 
spots. Longevity of star spots is at the heart of our approach. In 
order to obtain, in view of the extraordinary length of the time 
series, a satisfying fit, more freedom has been given to spot area 
evolution with respect to Paper 1 . It is now parameterized by up 
to eight parameters. 

Spot area is given in units of the star's cross-section. Area 
evolution is assumed to go basically linearly with time. The un- 
derlying physical reason is that then, at least in the case of a 
decaying spot, the slope of the area-time relation is somehow 
related to the turbulent magnetic diffusivity. With the aim of en- 
hancing flexibility and to describe the waxing and waning of a 
spot, three consecutive slope values are considered. The time 
derivative of spot area is then a mere step function over time. 
Step height measures the increase/decrease of area per day. So, 
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there are six free parameters: three slope values, two dates of 
slope change, and the logarithm of spot area at some point of 
the time series. We have done even a little bit more. In order to 
prevent sharp bends in spot area evolution, some smoothing is 
introduced. Each date where the slope changes is replaced by 
a time interval within which the slope is linearly interpolated 
between the two adjacent values. This makes the second time 
derivative of spot area a mere step function of time, described 
by six parameters. To get the integrated area itself as function of 
time two constants of integration enter, thus bringing the number 
of free parameters to describe a spot's area evolution to a total of 
eight. 

In addition to the free parameters of the model there are de- 
rived ones, the marginal distributions of which are of interest. 
An example is the rotational period. It follows from the longi- 
tudes of the spot centre at the beginning of the time series and at 
its end. 

All parameters are estimated in a Bayesian manner, i. e. their 
mean values as well as the corresponding uncertainties follow 
straightforwardly from the data alone. To maintain a flat prior 
distribution in parameter space, all dimensional parameters like 
periods or spot radii must actually be described by their loga- 
rithms. Only then the posterior probability distribution for a pe- 
riod will be consistent with that of a frequency and likewise the 
posterior for a radius with that of an area, i. e. it does not matter 
whether one prefers periods to frequencies or radii to areas . 

The likelihood function (Eq.O assumes that the measure- 
ment errors have a Gaussian distribution in the magnitude do- 
main. This is justified as long as the signal-to-noise ratio does 
not vary with changing magnitude, as it is for our data that span a 
full variation range of less than 0. 1 magnitudes. It has the invalu- 
able advantage that the likelihood function can be analytically 
integrated over measurement error cr, offset co, and linear trend 
do. To perfor m the integration over cr on e has to use Jeffreys' 
1 / cr-prior (cf . iKass & Wassermannl fl996h . The resulting mean 
likelihood depends on spot-modelling parameters pi ... Pm only. 
It takes into account all possible error values, offsets and lin- 
ear trends. By multiplying it with the prior, assumed constant in 
parameter space, one gets the posterior density distribution. All 
interesting quantities, parameter averages and confidence inter- 
vals, are then obtained by marginalization. 

With the N magnitude values d, measured at times f,, their 
standard deviations cr,, the model magnitudes /o(f;,/?i . . -Pm), 
offset cq, and trend do, the likelihood function is given by 



A(cr,CQ,do,pi ■ ..p M ;dd = 
(dj - f (VPum) - cp - dp ■ (tj - tp)) 2 

2cr? 



(3) 



1 J- ex p 

We set cr,- = s, ■ cr, with relative errors s,- being normalized ac- 
cording to l/sj =N. 

Parameter estimation by sampling the parameter space has 
bee n done by the Ma rkov chain Monte Carlo (MCMC) method 
(cf. lPress et al1l2007l) . 

Often parameter values are highly correlated. As MCMC 
performs best in an orthogonalized parameter space, all param- 
eters have been converted by a princi pal component an alysis us- 
ing singular value decomposition (cf. lPress et al.l2 007). Each pa- 
rameter in this abstract space is linearly dependent on all of the 
original parameters. The reconstruction of the original parameter 
values can be done exploiting a subspace of that orthogonalized 
parameter space. The dimension of that subspace, the number of 
degrees of freedom, proves lower by roughly one third or even 
more than the number of original parameters. 



4.3. Results 
4.3.1. KIC 7985370 

We have identified eleven gaps longer than an hour and two ad- 
ditional small jumps in the light curve (Fig. [8j. The data set was 
accordingly divided into 14 parts. Each part has been assigned 
its individual error level, off-set and, in Case B (i. e. non-rectified 
data), linear trend. Hence, the likelihood (Eq.|3) is the product of 
14 independent contributions. 

KIC7985370's inclination value i is - to be honest - ill- 
defined by the spot model applied to the Kepler photometry. 
Indeed, with only six spots the MCMC results in very dark spots 
(k « 0) at very low inclination (i « 10°). But even for these un- 
realistic solutions the equator-to-pole differential rotations was 
0.18 rad d~ 1 . Only with seven spots and allowing for enough spot 
evolution we arrived at acceptable inclination values (Pig.|9j and 
spot intensities (Fig. ITOb . If the inclination is fixed to the spec- 
troscopically derived value of i = 75° the residuals are rather 
high, ±2.46 mmag, exceeding the residuals of our best solution 
(+2.14 mmag) by far. Nevertheless, details of the solution with 
fixed inclination are also included in Table [5] where the results 
are presented. 

Improving the 7-spot solution by adding an eighth spot leads 
formally to a better fit. As the new spot proves to be ephemeral, 
lasting only six rotations, it neither constrains the differential ro- 
tation nor adds any significant insight (one can always get a bet- 
ter result by adding short-lived features). 



KIC7985370: inclination 




42 44 



Fig. 9. Determination of the stellar inclination from Kepler pho- 
tometry. Mean and 68-per-cent confidence region are marked 
by vertical lines (Case A only). Dashed: The corresponding 
marginal distribution for the original data with linear trends re- 
moved (Case B). 



The marginal distributions of the seven spot frequencies 
(Case A only), combined into one plot, are shown in Fig.fTTI 

From the three parameters describing the star's surface rota- 
tion, A, B and C, the equatorial rotational period (Fig. [12} and the 
equator-to-pole differential rotation (Fig. [T3l follow. The latter 
amounts to 0J774+°;°°{* radd -i (Case A) and 1729 + 0002 

radd 1 (Case B), respectively. The difference is significant, con- 
sidering the formal errors, albeit very small. In the case of fixed 
inclination (i = 75°) the differential rotation would be slightly 
enhanced, 0.1839 + 0.0002 radd 1 . 
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Fig. 8. Kepler light curve with best fit (solid red line, 2nd Case- A solution of Table [5} over-plotted. The residuals, shown at the 
top, are +2.14mmag. Obviously, the residuals are not homogeneous from one part of the light curve to another. 



KIC7985370: spot intensity 
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Fig. 10. Spot intensity related to the unspotted photosphere. 
Mean and 68-per-cent confidence region are marked by vertical 
lines (Case A only). Dashed: The corresponding marginal dis- 
tribution for the original data with linear trends removed (Case 
B). 



The spot area evolution is depicted in Fig. [14] The sudden 
rise of spot #7 seems to be an artefact. It falls into the gap be- 
tween the end of Q2 data and the beginning of Q3 data. On the 
other hand, the sudden disappearance of spot # 3 is not related to 
any switching from one part of the light curve to the next one. 
The fall in area of spot # 1 at the end of the time series is some- 
how mirrored in an increase in the size of spot #5. Maybe this 
indicates a flaw due to too much freedom in describing spot area 
evolution. 

Expectation values with l-cr confidence limits for various 
parameters are also quoted in Table[5] 

One should be aware that there is more than one solution 
for each case. The second Case-A solution presented in Table[5] 
is the one that has the lowest residuals found so far. There is 
an other well-relaxed 7-spot solution with slightly larger resid- 
uals nearby in parameter space. In that solution the fastest spot 
(#2), coming into existence near the end of the time series at 
JD~2455135, is located at a more southern latitude of -21°, re- 
sulting in a slightly increased differential rotation. All other spots 
are virtually unaffected. Further details of this second solution 
are given in Table[6] 



KIC7985370: spot frequencies 
5x1 o 4 t : 

4x10 4 r -. 

^ 3x10" r -. 
a. : : 

\ : : 

°- 2x10 4 r -: 
1 x 1 4 7 -. 

: - I . ; 

0.325 0.330 0.335 0.340 0.345 0.350 
1/P [d-] 

Fig. 11. All seven marginal distributions (Case A) of the spot 
frequencies. The three frequencies (0.324, 0.340, and 0.350 d" 1 ) 
seen in the low-resolution Fourier spectrum (Fig. |7) are con- 
firmed by the results of our spot model. 



Table 6. A second pair of 7-spot solutions for KIC 7985370. 



parameter 




Case A 


CaseB 


equ. period 


p 


2.8202 


±0.0002 


2.8209 


±0.0002 


diff. rotation 


dfi 


0.1943 


±0.0002 


0.1933 


±0.0002 


residuals 




±2.20 




±2.21 





Notes. The meaning of the entries is the same as in Table|5] i. e. peri- 
ods are in days, the differential rotation in radd~', and the residuals in 
mmag. 



4.3.2. KIC 77651 35 

We have identified eleven gaps longer than an hour and three ad- 
ditional small jumps in the light curve (Fig.[T5ll. The data set was 
accordingly divided into 15 parts. Each part has been assigned 
its individual error level, off-set and, in Case B (i. e. non-rectified 
data), linear trend. Hence, the likelihood function (Eq.[3j is the 
product of 15 independent contributions. 

As the inclination is photometrically ill-defined, we fixed it 
to the spectroscopically derived value of i = 75°. 
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Table 5. Three 7-spot solutions for KIC 7985370. Listed are expectation values and 1-cr confidence limits. 





parameter" 




Case A" 


Case A 1 ' 


Case B" 


inclination 


/ 


75!0 


fixed 


41?4 


+0.5 
-0.5 


32!6 


+0.7 
-o. 1 -) 


1st 


latitude 


R, 
Pi 


34!0 


+0.1 
-0.1 


29?4 


+0.4 
-0.4 


22!6 


+0.5 
-0.6 


2nd 


latitude 


o 
P2 


-10!0 


+0.6 
-0.9 


-6?9 


+0.9 
-0.9 


3?3 C 


+0.1 
-0.1 


3rd 


latitude 


O 
Pi 


32?2 


+0.2 
-0.2 


29?9 


+0.3 

-0.3 


21°TI 


+0,4 
-0.4 


4th 


latitude 


ft 
p4 


86!8 


+0.1 
-0.1 


87?5 


+0.1 
-0.1 


87.8 


+0.1 
-0.1 


5 th 


latitude 


a 

Ps 


53!8 


+0.1 
-0.1 


53?6 


+0.1 
-0.1 


51!1 


+0.2 
-0.2 


6th 


latitude 


P6 


35!8 


+0.2 
-0.2 


35?8 


+0.3 
-0.2 


33?0 


+0.3 
-0.3 


7th 


latitude 


a 
P7 


29!6 


+0.3 
-0.2 


19?9 


+0.9 
-0.9 


10!0 


+0.5 
-0.6 


1st 


period 


p 


2.8581 


+0.0001 
-0.0001 


2.8563 


+0.0001 
-0.0001 


2.8572 


+0.0001 
-0.0001 


2nd 


period 


p 


2.8350 


+0.0003 
-0.0003 


2.8428 


+0.0007 
-0.0006 


2.8475 


+0.0002 
-0.0002 


3rd 


period 


p 

^3 


2.8541 


+0.0004 
-0.0004 


2.8572 


+0.0004 
-0.0005 


2.8644 


+0.0007 
-0.0005 


4th 


period 


p, 
M 




+0.0001 
-0.0001 


J.U076 


+0.0001 
-0.0001 


J.Uoofi 


+0.0001 
-0.0001 


5th 


period 


Ps 


2.9382 


+0.0002 
-0.0002 


2.9421 


+0.0002 
-0.0002 


2.9417 


+0.0003 
-0.0002 


6th 


period 


p 

ih 


2.8629 


+0.0003 
-0.0004 


2.8700 


+0.0004 
-0.0003 


2.8754 


+0.0005 
-0.0004 


7th 


period 


P_ 
1 7 


2.8490 


+0.0003 
-0.0002 


2.8460 


+0.0002 
-0.0002 


2.8487 


+0.0002 
-0.0002 


spot intensity 


K 


0.437 


+0.005 
-0.004 


0.396 


+0.006 
-0.006 


0.406 


+0.012 
-0.012 


equ. period 


P 

1 eq 


2.8347 


+0.0003 
-0.0003 


2.8427 


+0.0007 
-0.0006 


2.8474 


+0.0003 
-0.0002 


deviation 


C/B 


0.28 




0.28 




0.21 




diff. rotation 


an 


0.1839 


+0.0002 
-0.0002 


0.1774 


+0.0004 
-0.0005 


0.1729 


+0.0002 
-0.0002 


residuals 




±2.46 




±2.14 




±2.12 





Notes. {a) Latitudes B are derived from the assumed law of differential rotation (Eq.[2}- Periods P are given in days, the spot intensity k is in units of 
the intensity of the unspotted surface. The ratio C/B measures the deviation from a pure sin~-law of differential rotation. The differential rotation 
dCl (radd~') is the equator-to-pole value of the shear. Residuals are in mmag. Case A refers to rectified data, Case B to non-rectified one. In 
order to get the Case-B solution the Case-A solution has been taken as a starting point for the MCMC parameter estimation. In the first Case-A 
solution the inclination is fixed to i = 75°. <c) The second spot is near the equator, therefore, the hemisphere it belongs to is doubtful. 



KIC7985370: equatorial period 
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KIC7985370: differential rotation 
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Fig. 12. Equatorial period of the star. Mean and 68-per-cent 
confidence region are marked by vertical lines (Case A only). 
Dashed: The corresponding marginal distribution for the origi- 
nal data with linear trends removed (Case B). 



Fig. 13. Equator-to-pole differential rotation of the star. Mean 
and 68-per-cent confidence region are marked by vertical lines 
(Case A only). Dashed: The corresponding marginal distribution 
for the original data with linear trends removed (Case B). 



Despite two spots more, the residuals, +2.35 mmag, exceed 
those of the seven-spot model of KIC 7985370 (+2. 14 mmag). 
This is not due to the fainter magnitude of KIC 7765 135 com- 
pared to KIC 7985370, because the photometric uncertainties are 



typically 0.047 mmag for the former and 0.022 mmag for the lat- 
ter. The reason may be that three of the nine spots are definitely 
short-lived with a life span as low as two months (cf. Fig. |20| >. 
which is less than twice the lapping time of 38 days between 
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KIC7985370: spot area evolution 



KIC7765135: spot frequencies 




1000 



1050 1100 
HJD - 2454000 
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Fig. 14. Spot area evolution (Case A). Area is in units of the 
star's cross-section. Vertical lines mark the boundaries of the Q0 
to Q3 quarters of data. A number in parenthesis indicates the 
spot number. 



the fastest and the slowest spot. We have to admit that dealing 
with nine spots goes already to the limit of the MCMC tech- 
nique since the method's relaxation time becomes prohibitively 
long. 

The marginal distribution of the spot rest intensity is shown 
in Fig. [16] 



KIC7765135: spot intensity 
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Fig. 16. Same as Fig.[TDJ for KIC 7765135. 



The marginal distributions of the nine spot frequencies (Case 
A only), combined into one plot, are shown in Fig.fTTI 

From the three parameters describing the star's surface rota- 
tion, A, B and C, the equatorial rotational period (Fig. [18} and the 
equator-to-pole differential rotation (Fig. [19} follows. The latter 
amounts to 0.1760 ± 0.0003 radd 1 (Case A) and 0.1774+°°^ 
radd" 1 (Case B), respectively. As for KIC 7985370 the differ- 
ence is small, but nevertheless significant. 

The level of differential rotation does not depend on the num- 
ber of spots considered. Neglecting the three short-lived spots, 
i.e. considering a six-spot model, would result in an equator-to- 
pole shear of 0.1777 + 0.0006 radd 1 . 



Fig. 17. Marginal distributions of the frequency for all the nine 
spots. The frequency values group around the two principal fre- 
quencies (0.391 and 0.414 d~') seen already in the Fourier spec- 
trum (Fig. [7}, which is the reason for the obvious "beating" phe- 
nomenon in Fig. [15] with a period of 40 days. The shortest and 
the longest frequency are a superposition of two frequencies. 

Inclination does not significantly affect dQ. Indeed, decreas- 
ing the inclination from the adopted value of i = 75° to 45° 
would result in a marginally larger differential rotation, three to 
four per cent. This is quite understandable. Inclination affects 
latitudes, but hardly periods. 

A cursory glance cast at the beating pattern (Fig. [15} re- 
veals a lapping time i\eat ~ 40 days, which is nearly exactly the 
lapping time of 40.3 days from the two peaks of the cleaned pe- 
riodogram (Fig. [7}. From these 40.3 days one already gets an 
estimate of the minimum value for the differential rotation as 
2^/Pbeat ~ 0.156radd~ 1 , which is not far from that derived by 
the model. 



KIC7765135: equatorial period 




2.4016 2.4018 2.4020 2.4022 2.4024 
P., [d] 

Fig. 18. Same as Fig. El for KIC7765135. 



The spot area evolution (Case A) is depicted in Fig. [20] The 
overwhelmingly large southern spot - at the beginning it fills 
to a large extent the southern hemisphere - may be an artefact. 
Because of its southern location its contribution to the light curve 
is rather modest. Perhaps it is actually a feature of the north- 
ern hemisphere, a non-circular extension of spot # 2. To prevent 



12 



H.-E. Frohlich et al,; Magnetic activity and differential rotation in two young suns 



KIC7765135 




950 1000 1050 1100 1150 

HJD - 2454000 

Fig. 15. Kepler light curve with best fit (solid red line, Case-A solution of Table [7]) over-plotted. The residuals, shown at the 
top, are +2.35 mmag. Obviously, the residuals are not homogeneous from one part of the light curve to another, hence, the stated 
+2.35 mmag is to be considered an overall average. 



KIC7765135: differential rotation 
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Fig.19. Same as Fig.Qj] for KIC 7765135. 

spot overlapping, spot # 8 had to be moved to the southern hemi- 
sphere. The reader should be aware that even in the case of a 
large spot the whole spot region has been assigned the angular 
velocity of its centre. Differential rotation is, to be exact, not 
compatible with a fixed circular shape. This is a shortcoming of 
our simple model. 

In the case of KIC 7765 135 it cannot be excluded that spot 
area evolution is partly driven by the need to avoid overlapping 
of spots. 

Expectation values with 1— cr confidence limits for various 
parameters are compiled in Table|7] 

5. Discussion 

5. ) . Chromospheric and coronal activity 

With the aim of making a comparison with the chromospheric 
activity of stars similar to our targets, we considered eigh- 
teen stars in the Pleiade s cluster, among those investigated by 
ISoderblom et al.l (1 199 3b). that have an early-G spectral type or 
an effective temperature close to our targets and a v sin i between 
4 and 18kms _1 . With a spectral subtraction analysis they found 
values of the net Ha emission between 100 and 500 mA and 
surface fluxes ranging from l.OxlO 6 to 2.9xl0 6 erg cm s . 



KIC7765135: spot area evolution 
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Fig.20. Same as Fig. [H] for KIC7765135. Three of the nine 
spots (#4, #7, and #9) are short-lived ones. Note the change in 
scale in the upper part! 



The net emission filling the core of the Ca H-/18542 line ranges 
from about 300 to 600 mA and the corresponding flux is in 
the range 1.0-2.7xl0 6 ergcm^s -1 . The values of these activ- 
ity indicators for our two targets are inside all these ranges, sug- 
gesting an activity level comparable to the Pleiades stars (age 
~ 130Myr). Th e X-ray luminosity fo r the se Pleiades m e mbers , 
as evaluated by IStauffer et all (11994 and iMarino et all (120031) . 
is in the range logLx = 28.7-29.7, which is just below the 
satura tion threshold of log Lx = 30.0 found by iPizzolato et al.l 
( 2003[) for Sun - like st ars. Moreover, according to th e works of 
Pizzolato et a D (120031) and ISoderblom etaD dl993bl) . the satu- 
ration would occur for rotation periods shorter than about 2 .0 
days. Another helpful work dealing with chromospheric activ- 
ity in stars belonging to the young open clusters IC2391 and 
IC2602 (age » 30 - 50Myr) shows that the chromospheric 
flux measured in the core of Can-/i8542 line saturates at a 
value of log(/?g 542 ) = logtFg^/oTeff 4 ) ~ -4.2 dMarsden et al.l 
I2009T). A similar behaviour was already found for the Pleiades 
bv ISoderblom et al.l (Il993bl) . The log(^^ 542 ) values of -4.5 and 
-4.6 that can be derived for KIC 7985370 and KIC 7765 135, re- 
spectively, are lower than, but not very far from, this saturation 
level. From the ROSAT X-ray count and the distance quoted 
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Table 7. Two 9-spot solutions for KIC 7765 135 with inclination 
being fixed to i = 75°. Listed are expectation values and 1-cr 
confidence limits. 



parameter 



Case A* 



Case B* 



inclination 


i 


75!0 


fixed 


75!0 


fixed 


1st 


latitude 


ft 


20?4 


+0. 1 

-d.i 


20°.0 


+0. 1 
-0.1 


2nd 


latitude 


ft 


76?3 


+(1.2 
-0.2 


75?1 


+0.2 
-0.2 


3rd 


latitude 


ft 


r.9 


+0.2 
-0.1 


0:8 


+0.2 
-0.2 


4th 


latitude 


Ih 


21 = 1 


+0.3 
-0.3 


20?5 


+0.3 
-0.3 


5th 


latitude 


ft 


0:9 


+0.3 
-0.4 


0:7 


+0.3 
-0.2 


6th 


latitude 


ft 


12?1 


+0. 1 
-0.1 


11?7 


+0. 1 
-0.1 


7th 


latitude 


Pi 


16?3 


+0.3 
-0.3 


17?2 


+0.3 
-0.3 


8th 


latitude 


ft 


-75!9 


+0.2 
-0.2 


-74!4 


+0.2 
-0.2 


9th 


latitude 


ft 


60!2 


+0.4 
-0.4 


59?4 


+0.4 
-0.4 


1st 


period 


Pi 


2.4223 


+0.0001 
-0.0001 


2.4222 


+0.0001 
-0.0001 


2nd 


period 


P 2 


2.5651 


+0.0001 
-0.0001 


2.5653 


+0.0001 
-0.0001 


3rd 


period 


Ps 


2.4020 


+0.0001 
-0.0001 


2.4022 


+0.0001 
-0.0001 


4th 


period 


Pa 


2.4237 


+0.0004 
-0.0004 


2.4231 


+0.0005 
-0.0005 


5th 


period 


Pi 


2.4019 


+0.0001 
-0.0001 


2.4022 


+0.0001 
-0.0001 


6th 


period 


A, 


2.4092 


+0.0001 
-0.0001 


2.4092 


+0.0001 
-0.0001 


7th 


period 


Pi 


2.4151 


+0.0004 
-0.0005 


2.4172 


+0.0005 
-0.0005 


8th 


period 


Ps 


2.5645 


+0.0001 
-0.0001 


2.5641 


+0.0001 
-0.0001 


9th 


period 


P 9 


2.5313 


+0.0010 
-0.0010 


2.5311 


+0.0009 
-0.0009 


spot intensity 


K 


0.671 


+0.003 
-0.003 


0.700 


+0.003 
-0.003 


equ. period 


P 

1 eq 


2.4018 


+0.0001 
-0.0001 


2.4022 


+0.0001 
-0.0001 


deviation 


C/B 


-0.008 




-0.008 




diff. rotation 


d£2 


0.1760 


+0.0003 
-0.0003 


0.1774 


+0.0003 
-0.0004 


residuals 




±2.35 




±2.29 





Notes. ' The meaning of the superscripts is the same as in Table [5] 



in Table|2] we evaluated the X-ra y luminosity of KIC 7 985370 
through the relation proposed by iFleming et al.l (1 19951) . Lx = 
And 2 {%3\ + 5.30 H R1) x 10~ 12 ergs"' , where the hardness ra- 
tio is HR\ = 0.23 (IVoges et al.ll2000l) . The value of logL x = 
29.68+° ^ confirms the non-saturated regime for KIC 7985370. 



5.2. General considerations on the spot model 

Despite the fact that both stars are very active ones and exhibit 
filled-in absorption of several chromospheric activity indicators, 
our photometric analysis is in terms of dark surface features only. 
Allowing for bright ones too, would make the MCMC approach 
usually unstable. Anyway, photometry alone seems to be un- 
able to discriminate even between dark and bright spots (e.g., 
lLuftingeretaT1l2010l) . 

As a comparison of the two cases A and B reveals, the results 
of our spot modelling are hardly influenced by the rectification 
procedure. The smaller residuals for non-rectified data (Case B) 
are very likely due to the fact that the Case-B likelihood func- 
tion (Eq.0 also takes into account a linear trend in the data, in- 
dividually for each part of the light curve. This allows for more 
freedom in fitting the data and results in a slightly better fit. 

The reader should be aware that the estimated parameter val- 
ues and their often surprisingly small errors are those of the 



model constrained by the data. Error bars indicate the "elbow 
room" of the model, nothing more. 

5.3. Frequencies 

It is remarkable that the frequencies that stand out in the power 
spectrum of the light curves (Fig. [7} represent the distribution 
of spot frequencies (Figs. QT| and ITTb astonishingly well. The 
lapping time, as a measure of the lower limit of surface differ- 
ential rotation, follows already from the periodogram analysis! 
However, in order to get an estimate of the full equator-to-pole 
span of the latitudinal shear, including its sign, one needs latitu- 
dinal information. 

5.4. Inclination 

Combining the inclination value from photometry, i * 40°, 
with the spectroscopically measured projected rotational ve- 
locity vsin/ (Table |2]i allows us to determine in the case of 
KIC 7985370 the star's radius. Taking the shortest rotational 
period (P eq ), one arrives at R = 1.42-1.44J?q. This (minimal) 
radius is larger than the ZAMS value of R w 1.1 Rq, but smaller 
than the R * 2Rq for a star of 1.5M© at 10 Myr. Hence, the 
photometrically derived radius is within the expected range. As 
stated in Sect. 14.3.21 in the case of KIC 7765 135 the photometric 
inclination is badly defined. Therefore, the inclination has been 
fixed to i = 75°, assuming the radius to have its ZAMS value. 

5.5. Spot contrast and spot longevity 

Although both stars share the same spectral type and age, there 
are differences concerning the spots. The spots of KIC 7985370 
seem to be darker and longer living than those of KIC 7765 1 35 . 

We would like to remind that a "spot" may be in fact a group 
of smaller spots that all together form an active region, which 
could also include bright features. 

For KIC 7765 135, the spot contrast, k w 0.7, looks 
rather normal. It is similar to the previously studied case of 
KIC 8429280 (Paper I). The corresponding temperature contrast, 
the ratio between spot and photospheric temperature r sp /r p h, 
is 0.9, assuming that the "white light" Kepler flux matches the 
bolometric conditions. The much darker spots, k » 0.4, in the 
case of KIC 7985370 defy a simple explanation. There is no need 
for exceptionally small (and therefore dark) spots to prevent spot 
overlap. 

Apart from a few late F-type stars with lo w or moderate ac- 
tivity observed by CoRoT dMosser et al.l2009l) where spots seem 
to be short-lived, there is strong evidence that spots in very active 
stars like our targets have rather long lives compared to the star 
rotation. Active longitudes lastin g for months or years have been 



observed in young stars (e.g., Collier Cameronl 1995 ; lHatze"sl 
Il99l lBlrnes^Lal][T99l Huber et al.ll2009t lLanza et al.l 1201 ll) 
and in the evolved componen ts of close bin ary sy stems, like 
II Peg (e.g. lRodono et al-lhoOCtlLindborg et all201 ll) . This does 
not exclude that individual unresolved spots, which are compos- 
ing the active region, have shorter evolution times, but the pho- 
tospheric active region seen as an entity endures for a very long 
time in such cases. 

Unlike KIC 7985370, in the case of KIC 7765 135 mid- 
latitude spots (30-50°) are missing. This is reminiscent of the 
spot distribution of two fast-rotating early G dw arfs, He 520 
and H e 699, of the a Persei cluster studied by iBarnes et al.l 
(Il998h . Despite the fact that there is a clear distinction between 
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near-equator and near-pole spots, with regard to spot lifetimes, 
no correlation seems to exist between lifetime and latitude, 
whic h is contrary to the cas e of the rapidly-rotating young AB 
Dor dColher Cameronfl995l) . where only low- and intermediate- 
latitude spots are long-lived. 

The sudden appearance of a full-grown spot at the beginning 
of a new quarter of data is suspicious and hints at an artefact 

(Fig.na. 

5.6. Differential rotation 

Both stars exhibit low-latitude spots as well as high-latitude ones 
at the time of observation making them suitable for studying 
their latitudinal shear. The most robust and important result of 
the present work is the high degree of surface differential rota- 
tion found for both stars: dQ = 0.18radd -1 . This exceeds three- 
fold the solar value. 

This estimate is rather robust, because any spot model with a 
few long-lasting spots able to reproduce the beating of the light 
curve must provide a value of equator-to-pole differential rota- 
tion that exceeds the lower limit of irrespective of the 
number of spots used. 

Inclination has a marginal effect since the periods found in 
the light curve do not depend on it. 

We remark that the high value dQ relies on the assumption 
of spot longevity. It is always possible to get an excellent fit with 
many short-lived spots even for rigid rotation. 

Very different values of differential rotation have been found 
for HD 171488 (V889 Her), a young (-50 Myr) Sun. For this 
star, which is rotating faster (P = 1.33 days) than our tar- 
gets, a very high solar-type differential rotation dQ ss 0.4-0.5 
radd -1 , with the e quator lapping t he po les e very 12-16 days, 
was fo und by both Mar sden et al.l (120061) and iJeffers & Donatil 
(12008b . Much weaker v alues (dQ » 0.04) were derived in- 
stead f o r the same st a r by prvinen etall (120081) and lKovari et alj 
(1201 ll) . iHuber et al.l i2QQ% even claim their data being consis- 
ten t with no differential r otation. 

iMarsden et al.l d201 lb report values of dQ in the range 0.08- 
0.45 radd -1 for a sample of stars similar to and slightly more 
massive than the Sun. Among these stars, HD 141943, a 1.3- 
Mq star that is still in the PMS phase (age ~ 17 Myr), displays 
values of dQ ranging from about 0.23 to 0.44 radd -1 in differ- 
ent epochs. A sola r-type differential r otation, dQ « 0.2 radd -1 , 
was also found bv lWaite et all d20Tlb for HD 106506, a Gl V- 
type star (r e ff = 5900 K) that is very similar to our targets, but it 
is rotating faster (P e q = 1.39 days). Moreover, the Fourier trans- 
form technique applied to high-resolution spectra of a large sam- 
ple of F- and early G-type sta rs indicates that differen t ial rota- 
tion i s rather frequently found dReiners & Schmittl2003l iReinersI 
[2006b . In their data, there is no clear dependency on the rotation 
period, but the strongest differential rotation, up to ~ 1.0 radd -1 , 
occurs for periods between 2 and 3 days and values as high as 
~ 0.7 radd -1 are encountered down to P ~ 0.5 days. 

From ground-based photometry, which is basically devoted 
to cooler stars, a different behaviour, i.e. a differential rota- 
tion decreasing with the rotation period, seems to emerge (e.g., 
M essina & Guinanl2003l) . However, the precision of the ground- 
based light curves does not allow to draw firm conclusions and 
accurate photometry from space, as well as Doppler imaging, is 
needed for settling this point. 

For mid-G to M dwarfs, weaker v alues of the latitud inal 
shears are generally found. In particular. iBarnes et al.l(l2005l) an- 
alyzed with the Doppler imaging technique a small sample of 
active stars in the spectral range G2-M2 finding a trend towards 



decreasing surface differential rotation with decreasing temper- 
ature. This suggests that the stellar mass must also play a sig- 
nificant role in this respect. The largest values for stars as cool 
as about 5000 K are dQ = 0.27 radd -1 found by us in Paper I 
for KIC 8429280 (K2 V, P = 1.16days) and dQ = 0.20 radd" 1 
found bv lDonati et all J20031) for LQ Hya (K2 V, P = 1 .60 days). 
The slowly rotating (P e q = 11 -2d) and mildly active K2V star 
eEri exhibits only little differential surface rotation (0.017 < 
dQ < 0.056 radd -1 ) as a Bayesian r eanalysis of the MOST lig ht 
curve (ICrollll2006HCroll et al.ll2006b revealed (lFrohlichl2007b . 

Thus, there is an indication that a high differential rotation 
goes along with a high rotation rate. 

The high differential rotation that we found for KIC 7985370 
and KIC 7765 135 disagrees with the hydrodynamical model of 
iKiiker et"al] d201 ll) . which instead predicts a rather low value of 
dQ w 0.08 for an (evolved) solar-mass star rotating with a period 
as short as 1.3 days. 

Surface differential rotation may even vary along the activ- 
ity cycle. Indeed, certain mean-field dynamo models for rapidly 
rotating cool stars with deep convection zones predict torsional 
oscilla tions with variation s of several percent in differential ro- 
tation dCovas et al ] l2005h . Of course, this cannot explain such 
extreme cases as LQ Hya where at time s the surface rotation is 
solid body. According to lLanz al d2"006h . to maintain the strong 
shear (~0.2 radd -1 ) observed for LQHya in the year 2000 
would imply a dissipated power exceeding the star's luminosity. 

The differential rotation of rapidly-rotating solar-like stars 
has been recently inves tigated on theoretical grounds by 
iHotta & Yokovama1d201 lb . They found that differential rotation 
approaches the Taylor-Proudman state, i.e. the iso-rotation sur- 
faces tend to become cylinders parallel to the rotation axis, when 
stellar rotation is faster than the solar one. In this case, the dif- 
ferential rotation is concentrated at relatively low latitudes with 
large stellar angular velocity. They show that the latitudinal shear 
(between the equator and latitude /3 = 45°) increases with the an- 
gular velocity, in line with our results and the recent literature. 



6. Conclusions 

In this paper, we have studied two Sun-like stars, KIC 7985370 
and KIC 7765 1 35, by means of high-resolution spectroscopy and 
high-precision Kepler photometry. 

The high-resolution spectra allowed us to derive, for the first 
time, their spectral type, astrophysical parameters (T e g, \ogg, 
[Fe/H]), rotational and heliocentric radial velocities, and lithium 
abundance. All this information, combined with the analysis of 
the SED and proper motions, enabled us to infer their distance 
and kinematics, and to estimate the age of both stars in the range 
100-200 Myr, although we cannot exclude that they could be as 
young as 50 Myr. Thus, these two sources should be already in 
the post-T Tauri phase. 

As expected from their young age, both stars were found 
to be chromospherically active displaying filled-in Ha, H/3, and 
Ca n IRT lines, as well as He i D3 absorption. The surface chro- 
mospheric fluxes and the X-ray luminosity (for KIC 7985370), 
within the ranges found for stars with similar r e ff and vsinz 
in the Pleiades cluster, are just below the saturation level 
dSoderblom et al.lll993bl) . The flux ratio of two Ca 11 IRT lines 
and the Balmer decrement (for KIC 7765 135 only) suggest that 
the chromospheric emission is mainly due to optically-thick sur- 
face regions analogous to solar plages. 

We have applied a robust spot model, based on a Bayesian 
approach and a MCMC method, to the Kepler light curves 



15 



H.-E. Frohlich et al,; Magnetic activity and differential rotation in two young suns 



which span nearly 229 days and have an unprecedent preci- 
sion (* 1(T 5 mag). While seven long-lived spots were needed 
to perform a reasonable fit (at a 2-mmag level) of the light 
curve of KIC 7985370, we used up to nine spots in the case 
of KIC 7765 135 due to a shorter lifetime of its spots. Because 
of the exceptional precision of the Kepler photometry it is im- 
possi ble to reach the Bayesian noise floor defined by, e. g., the 
BIC dSchwarz|[l978T) without increasing significantly the de- 
grees of freedom and, consequently, the non-uniqueness of the 
solution. Provided spots are indeed long-lived, the equator-to- 
pole value of the shear amounts for both stars to 0.18 radd 1 . 
This i s in contrast with the theoretical models of iKiiker et alj 
(1201 ll) that predict a moderate solar-type differential rotation 
even for fast-rotating main-sequence stars, unless the convec- 
tion zone is shallower than predicted by the stellar models. Our 
results are instead in line with the scenario proposed by other 
modelers of a differential rotation that increases with the angu- 
lar velocity dHotta & Yokovamall201 ll) and that can be also sub- 
ject to changes along the activity cycle (ICovas et alj2005HLanzal 
120061) . 
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